function [ W1 H1 ] = GetStartingWHMySVD( duomenys, R )
    [U,S,V] = svds(duomenys);
    
    W1 = U * S;
    H1 = S * V';
    
    W1 = W1(:, 1:R);
    H1 = H1(1:R,:);
    
    ee = 0.1;
    
    W1(W1 <= 0) = ee;
    H1(H1 <= 0) = ee;
    
    W1 = W1 / max(max(W1));
    H1 = H1 / max(max(H1));
    
end

